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The general principles needed to compute the effect of a stationary gravitational field on the quasistationary elec- 
tromagnetic phenomena in normal conductors and superconductors are formulated from general relativistic point of 
view. Generalization of the skin effect, that is the general relativistic modification of the penetration depth (of the 
time-dependent magnetic field in the conductor) due to its relativistic coupling to the gravitational field is obtained. 
The effect of the gravitational field on the penetration and coherence depths in superconductors is also studied. As 
an illustration of the foregoing general results, we discuss their application to superconducting systems in the outer 
core of neutron stars. The relevance of these effects to electrodynamics of magnetized neutron stars has been shown. 
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1. Introduction 

Since electromagnetic system is a purely relativistic system, its coupling to gravity can be described, in prin- 
ciple, only general-relativistically. In our previous research, the effect of the gravitational field on stationary 
electromagnetic processes in conductors and superconductors in the presence of constant electric current has 
been investigated. In fact, the influence of angular momentum of the gravitational source may appear as a 
galvanogravitomagnetic effect in current carrying conductors, 1 as general-relativistic effect of charge distri- 
bution inside conductors in an applied magnetic field, 2 as general-relativistic effects in charged systems. 3 
The effect of gravitational field on thermoelectric effects in superconductors and on electromagnetic prop- 
erties of type II superconductors have been investigated in our papers. 4 ' 5 Here we extend our results to the 
time-dependent electromagnetic processes. 

The paper is organized as follows. First we write Maxwell equations in the Schwarzschild space-time in 
Section 2. Then in the next part of the paper, we will study the effect of gravity on the quasistationary 
time-varying electric current and magnetic field penetrating the conductor placed in the gravitational field. 
Sections 3 is devoted to investigation of the skin effect in the Schwarzschild space-time and in a weak gravi- 
tational field when the curvature effects are locally negligible. In particular, new purely general-relativistic 
effects connected with the gravitational modification of penetration depth of magnetic field will be obtained. 

The remarkable quantum phenomena exhibited by a superconductor show that the superconducting state 
in Ginzburg-Landau theory 7 described by a wavefunction can be spread over the entire superconductor. The 
existence of macroscopic wavefunction raises the natural question of how it is influenced by the gravitational 
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field. In the last part of the paper, that is, in Sections 4-7 we shall formulate the general formulae and their 
solutions needed to determine the influence of stationary gravitational field on the quasistationary electro- 
magnetic phenomena for superconductors, from general- relativistic point of view and study the relevance 
of general-relativistic effects in superconductors of II type on electrodynamics of magnetized neutron stars. 
This analysis has been motivated by the mixture commonly found in neutron star models, namely in the 
outer core region, where supcrfluid neutrons, superconducting protons and normal electrons are generally 
thought to coexist (see for example, Ref. 8). However, due to the generality of the present approach, it is 
equally well applicable to superconducting systems found in more common laboratory contexts in the weak 
gravitational field of Earth. 

In this paper we use a space-like signature (—,+,+,+) and a system of units in which G = 1 = c 
(however, for those expressions with a physical application we have written the speed of light explicitly). 
Greek indices are taken to run from to 3 and Latin indices from 1 to 3; covariant derivatives are denoted 
with a semi-colon and partial derivatives with a comma. 



2. Maxwell Equations in Schwarzschild Space-time 

In a coordinate system (t, r, 9, <p), the metric for a static spherically-symmetric relativistic star with the total 
mass M is 

ds 2 = -N 2 dt 2 + N~ 2 dr 2 + r 2 d9 2 + r 2 sin 2 9d<p 2 , (1) 

where N 2 = 1 — 2M/r is the lapse function. 

The general form of the first pair of general relativistic Maxwell equations is given by 

Fa/3, 7 + Fja,0 + Ff3y,ct = , (2) 

where F a p is the electromagnetic field tensor expressing the strict connection between the electric E a and 
magnetic B a fields. For an observer with four-velocity (u a ) obs , the covariant components of the electromag- 
netic tensor are given by 

F aSi = 2(u) obs[a E m + r, a ^ 5 (u^ obs B s , (3) 
where A[ a ^ = ^(A a p — Ap a ) and rjap^s is the pseudo-tensorial expression for the Levi-Civita symbol e Q/ 3 7< 5 

l] af3jS = - -j== e a 0jS , Va/S^S = \/~-gtaSi 1 8 , (4) 

with g = det\g a p\ — — r 4 sin 2 9 for the metric (1). 

The proper observers have four-velocity components given by 

(u a ) piop = N- 1 (l, 0, 0, o) ; K) prop = N ( - 1, 0, 0, o) . (5) 

In the metric (1) and with the definition (3) referred to the observers (5), the first pair of Maxwell 
equations (2) take then the form 11 

sin9(r 2 fT) r + £ (stadB^ + = , (6) 
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(^ sm ")^ = sm9 ( rW£ *), r - E '*< (8 » 

-(""%■ < 9 > 

where "hatted" quantities are projected onto a locally orthonormal tetrad {e^} = (e^, e r , e^, e^) carried by 
a proper observer 

e? = l(l,0,0,o) , e? = Jv(o,l,0,o) , ef = \ (0, 0, 1, 0) , e| = ^ (0, 0,0, l) . (10) 

The general form of the second pair of Maxwell equations is given by 

F Q/3 . p = 4tt J a (11) 

where the four-current J a is a sum of convection and conduction j a currents 

J a = Pe u a +j a , j a u a = 0, (12) 

with u a being the conductor four-velocity as whole and p e the proper charge density. 
We can now rewrite the second pair of Maxwell equations as 11 

(s,„ob*) s -s%= (Tfl^ + i^er , (i4) 

B'\ t - sin 6 (rWB*) = (" ! ^) ^ + 4»r sin 0.? , (15) 
3. Skin Effect in Schwarzschild Gravitational Field 

If the conduction current j a is carried by the electrons with electrical conductivity a, Ohm's law can then 
be written as 

j a = aF^u 13 . (17) 

Equations (14)-(16) can now be rewritten in a more convenient form. Taking four-velocity components 
of the conductor as (5), we can use Ohm's law (17) to derive the following explicit components of current 
J a in the proper frame 

J* = Pe , ,f — oE l . (18) 

Hereafter we will use the following realistic assumptions. Firstly, we consider a to be uniform within the 
conducting media. Secondly we ignore the contributions coming from the displacement currents. The latter 
could, in principle, be relevant in the evolution of the electromagnetic fields, but their effects are negligible 
on timescales that are long as compared with the electromagnetic waves crossing time. In view of this, we 
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will neglect in (14)-(16) all terms involving time derivatives of the electric field (quasistationarity condition) 
and use Ohm's law (18) to rewrite Maxwell equations (14) and (15) as 



EP 
E 4> 



Airar sin 9 

c 

Anar sin 9 



sm 



B f t4> - sin e(rNB^ 



(rNB 8 ^J - B f 



(19) 
(20) 
(21) 



Using Maxwell equation (13) and Ohm's law (18), we find that the space charge density p e = p e (t, r, 9, <j>) 
inside the star has a zeroth-order contribution given by 



Pe = 



cN 



Airr 2 sin 9 



sin 6{r 2 E r ) r + ^ (sin «/•;*) - — E' 



r 

,o ' N' 



(22) 



Now with the help of Eqs. (7)-(9), (19)-(21) and Ohm's laws (18), we obtain the induction equations for 
the components of time varying magnetic field penetrating the conductor 



dB r 
~df 

dB § 
dB^ 



c z N 



1 



47rcrr 2 sin 9 I sin 9 
c 2 N 



B' 



Anar sin 9 



sin 6 N 



<p -sm9(rNB< 
(rNB §S ) - B\ 



— sm( 



r sm( 



c 2 N j 

47rcrr 1 \ r sin 



B« 



N 



sint 



B\ 



(rNB°) r -B^ 
(smeB^ ^ - B § 
- sin 9 (rNB^j 



(23) 
(24) 
(25) 



Assume a conducting sphere with the conductivity a occupies the space r < R with empty vacuum space 
in the region r > R. The surface at r = R is subjected to dipolar magnetic field varying harmonically with 
time. For example, the azimuthal component of magnetic field B e is 



B e = B^eT 1 " 1 



(26) 



In order to understand the behavior of the magnetic field inside the spherical conductor we look for a 
solution of the induction equation (24) in the inner region r < R 

c 2 N 



-iuB° 



Airar 



N 



rNB 



which can be written as the following differential expression 



d 2 B B 
dr 2 



3M 

r 2 ]\[2 



2\ dB e 
dr 



2 130 



k z B 



M 



B» = 0, 



(27) 



(28) 



where the parameter k 2 = (Airaiui) / (c 2 N 3 ) is introduced. 

Analytical solution of this equation can not be found easily, and as we do not need any numerical solution 
for our purpose we must choose some approximation that do not destroy the structure of possible solution. 
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Due to this purpose we will see that in linear M/r approximation this equation could be simplified and 
written as Bessel equation 



d 2 B° . 2dB* . , 2 _<| 



+ k?B" = . 



(29) 



dr 2 r dr 

The boundary condition for this equation is B® — B e (R). And the solution of this equation is the spherical 
Bessel function of the first type 12 

,§ Jo(fcr) 



B" = b: 



°jo(kR) 



where 



k = 



1 + i 



5 = 



cN 3 ' 2 



(30) 



(31) 



o \f2iujuj 

We will expand the expression for the magnetic field (30) at the surface of the sphere for the large 
arguments of kr » 1. If we take i = i?-rasa distance from the surface inwards then the interior magnetic 
field behaves as 



B e = B e e ikx = B e e if + B o e -f ; 
and the distance x — S is called the skin depth of magnetic field penetration (see Fig. 1). 



(32) 




0.002 



0.004 O.OOfi 

z a 



o.oos 



0.01 



Fig. 1. The dependence of skin depth in the proton stellar superconductor on the compactness parameter i n the case 
when conductivity a »J 10 25 s _1 and the frequency ui fs 10 3 s _1 . 



From Eq. (21) one can find the expression 



i _ ckN ji (kr) § 
47ra j (kR) *° ' 



(33) 



for azimuthal component of the electric field. 

Thus the electric field falls off exponentially in r, with a spatial oscillation of the same scale, being 
confined mainly to a depth less than the relativistically modified skin depth. 
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3.1. Skin effect in weak Earth's gravitational field 

In the weak gravitational field of the Earth the curvature effects are negligible and space-time metric in the 
Cartezian coordinates is 



ds z = -c 2 N 2 dt 2 + dx 2 + dy 2 + N-' z dz z . 



2j,2 



Here N is equal to 



N = 1 



2gz 

~~2~ 



(34) 



(35) 



and z is the height above some fixed point, when the apparatus at the Earth may be regarded as having an 
acceleration g relative to local inertial frame. 

The evolution equations for the magnetic field in the metric (34) take the form 



dB £ 
dt 


Aira 


dBy 

dt 


c 2 

= N 

Ana 


dB z 
dt 


= -^—N 

Ana 



N 



(NB & ) 



B z 



N 



Bf v - (.v/r'S 



-[ B %- B %], y 

Bf y - (NB V ) 
(NB*) -B% 



(36) 
(37) 
(38) 



Assume a semi-infinite conductor of uniform conductivity a occupying the space z > 0. The surface at 
z = is subjected to tangential magnetic field B x = B^e" 1 ^ governed by the differential equation 

(39) 



d 2 B* 6g dB* Ag 2 
-^ + c^N^z- + c^P B +kB -°' 



which has a solution 

_ ( j^ e (-3c 2 gN-y / 5c 4 g 2 N 2 -c 8 k 2 N 4 -)z/c 4 N 2 _^ (7 2e (~ 3c2 9 N +V 5c 4 fir 2 N 2 -c s k 2 N^z/c 4 N 2 j-^q-j 

where C\ and Ci are the integration constants. Neglecting small terms and using (32) for k one can get 



4. Effect of Gravity on Coherence Length and Penetration Depth 

The coherence length is the measure of the distance within which the properties of the superconductors are 
not changed appreciably in the presence of a magnetic field. Because we are concerned with the effect of 
gravity on the coherence length in a superconductor it would be useful first to discuss Ginzburg-Landau 
equation in a gravitational field 7 



which could be written as 



aip + (3ip\ip\' 



(42) 
(43) 
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Here tp is the complex order parameter describing the superconductor, a and (3 are the phenomenological 
expansion coefficients. 

We now rewrite Ginzburg-Landau equation (42) in the Schwarzschild metric (1) in a more useful form 



1 f 



4m 
+2ei 

+2ei 
+4e 2 



JPdt 2 

1 dAt 



dr 



r 2 N 2 - 



d_ 

09 



sin 9- 



1 d 2 



N 2 dt 

A! 

N 2 dt~ 



+ 



J_ (sinO^- [r 2 N 2 A r ] + A [ s i n 0A e ] + 
smO \ or 1 J aO 



sin# 
1 <9A 



sin# d(j) 



i 



i 



4 2 



>V = o . 



(44) 



Assuming that the order parameter ip has no variation in coordinates 9 and <j> and all components of 
vector-potential of electromagnetic field vanish then Eq. (44) reduces to 



1 



«V> + VMM" ~ -^N 2 ( 



(d 2 i> 



\ dr 2 ^ r 



1 + 



M 

Vn 2 



= o 



which in approximation Mj r«l takes the simple form 



4m 



2 -JLjV 2 {^ + ^± 



^ dr 2 ^ r dr 



= 



Now we introduce the dimensionless wavefunction ip(r) 



<P(r) = 



where 



2 



\a\ 



C 2 (^ + 2 -^)=o 

\ dr 2 r dr J 



(45) 
(46) 

(47) 

(48) 
(49) 



Here the parameter £ = £, 2 N = 4 ^\a\ ' £° ^ s ^ ne Ginzburg-Landau coherence length in the flat Minkowski 
space-time, n( s ) is the density of superconducting electrons. The coherence length is the length scale on which 
the condensate charge density 2e|-0| 2 vanishes, as one approaches the surface of the superconductor from its 
interior. 

Now we distinguish the meaning of £. The coordinate axis r is perpendicular to the surface of the 
superconductor r = R. Assume that the surface layer is slim and the value of the function tp at the surface 
has small deflection from 1, that is 



<p = l- e(r) , 



(50) 
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and substituting the expression for ip into (49) we obtain the Bessel equation 

— + -— -2C 2 e-0 
dr 2 r dr 



12 



(51) 



in the linear approximation in e. 

Then taking into account the property lim r ^ <P = 1 one can get e(0) = and find the solution of the 
Eq. (51) as modified Bessel function of the first type 12 



e = e(0) — 



jo(iv^f) 



jo(iv^f) ' 



(52) 



It is known that the spherical Bessel function j (r) of zeroth order for large arguments r » 1 is approx- 
imated as 

sinr 



jo(r) 

i 

and due to this approximation near the surface of the sphere expression (52) takes the form 

e = e(0)e -V5 f , 



(53) 



(54) 



where x — R — r is the inward distance from the surface of the sphere. It follows that £ is the length scale 
of the parameter of order <p (See Fig. 2). 




4X10 



3X10 



0.2 



0.4 



0.6 



0.S 



Fig. 2. Exponential growth of the order parameter with the distance x = R — r from the surface of the gravitational object 
toward its interior; normal line is responsible for Newtonian fiat space-time and dashed line is for neutron star with the 
compactness parameter 2M/r = 0.5. 



In the relativistically generalized London theory the equation for the superconducting current jf is 6 



2j( s )[ Q; /3] - 2j( s )[ a (lnn( s )) )(8 ] + cp ois )(Ap a + 2u[ a w ]) 



4n 



-XV'F, 



a/3 



(55) 



where the notation A| = (mc 2 ) / (Ann^e 2 ) for the London penetration depth is introduced and 

A a p = u [a , fj] + u [fj w a] (56) 
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is the relativistic rate of rotation, w a = Ua^u 13 is the absolute acceleration of superconductor. 

Using the London equation (55) we obtain the expressions for the electric field (we choose c=l for this 
case too) 

E f = AirXllN- 1 ^ - (Nf), r - AT 1 / (Inn),* + JV/(Inn), r + P N r ] , (57) 



E e = 4 7 rA|[AT- 1 4 - r- 1 ^ - JV" 1 / (lnn), t + r-V(Inn), fl ] , 



(58) 



E* = 4 7 rA|[Ar- 1 4 - r-^in" 1 ^ - N~ x f (lnn), t + r" 1 sin" 1 0/(1™),^] . (59) 

Here the derivatives from the density n( s ) of superconducting pairs are taken to be equal to zero, because 
the dependence of ri( s ) on the coordinates is comparably small. Putting f extracted from the second pair 
of the Maxwell equations (14)-(16) into expressions (57) (59) in the framework of this approximation and 
inserting the obtained result into the first pair of Maxwell equations (6)-(9) one obtains the expressions for 
magnetic field which are the second London equations 



1 



rsmv rsinf 



X\N 



rsin( 



B^-sin( 



rNB 



sm( 



B\ 



,t<t> 



1 



rNB 



B\ 



Nr sin 9 



sin 



9B 4, 



B ti 



1 



Nr sin 9 



sin(9B^ 



-B. 



1 



AO 



sinf 



,t<t> 



-NB^ 



For simplicity of calculations the London equation (61) for B e component of the magnetic field 

B> - (rNB') . A|i^ (-^ + H B » + + flj 

can be approximately written as Bessel like equation 

2 



rr T D r — A £> 

r ' 







in the linear approximation in the small compactness parameter M/r, where parameter A = X^N. 
The solution of the Eq. (64) is already known and expressed through the Bessel function 



B ti = B[ 



9 Jo (» I ) 
%(*f) 



(60) 

(61) 
(62) 

(63) 
(64) 

(65) 



Using the above mentioned arguments one may get the following expression at the surface of the sphere 
(see Fig. 3) 



B° 



(66) 



Our main finding is the dependence of the coherence length and penetration depth on the gravitational 
field. It is connected with the general relativistic modification of surface electromagnetic energy of the 
conductors and superconductors due to the effect of gravitational field. 
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Fig. 3. The exponential decay of magnetic strength inside a superconductor of type II, normal line is for flat space-time, dash 
line is for a neutron star with the compactness parameter = 0.5. Here x = R — r is the distance from the surface of the 
star toward its interior. 



5. Magnetic Field of Single Vortex 

The four-vector of supercurrent density is 



3{s)a — 



le 



iP*(d a - 2eiA a )i) - 4,{d a + 2etA a )i>* \ , 



(67) 



which could be written through dimensionless function (p — ip/ipo as 



3(s)c 



- 1 — — ((fi*d a tp - fd a ip*) - —^A r _ 



(4ttA) 2 



4ttA 2 ' 



where 



he 
2e 



is the quantum of magnetic flux and 



A 



m(3c 2 



(68) 



(69) 



87re 2 |a| 

Assume that the wavefunction is ip — \ f\e li> . Then superconducting current density (68) takes the form 



(70) 



3(s)c 



4ttA 2 V2tt 



<I>n 



1 f ^-A a 



4ttA 2 V 2tt 



(71) 



If the parameter of the Ginzburg-Landau theory n — A/£ is much bigger than 1 then at large distance r>( 
we will have \ip\ 2 ~ 1. 

In cylindrical coordinates (t, r, <fi, z) the space-time metric around a gravitating source takes the form 10 



ds 2 



N 2 dt 2 +N- 2 (dr 2 



dz 2 



(72) 
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Using this metric we take curl from both sides of Eq. (71) and obtain the differential equations for the 
components of magnetic field B l : 



b: 



B: 



rB" 



According to Stokes Theorem 13 



b: 



[ B l'zz - B *rz] 

\rB^). r -B% 
rB* 



X 2 N 2 | 2tt 
1 



(V x Wf - B f 



— (V x - -B* 



X 2 N 2 2ir 



X 2 N 2 2ir 



(V x w) £ - B z 



(V x W)" = 2ir5(x a - afljfj" , 



(73) 
(74) 
(75) 

(76) 



here rf is the unit vector aligned along the vortex. 

Assume B^ = B r = and B z {r) ^ 0, then Eq. (75) could be written as equation for the z component of 
magnetic field 



B z - X 2 N 2 | B 



1 



-B 2 



(77) 



where e z is the unit vector along z axis. This is the Bessel equation of complex argument with the boundary 
condition B z (oo) = 0. The solution of the Eq. (77) is the Bessel function of the third type, either the 
McDonald function or Hankel function of complex argument 

$ 



B z 



;K (r/XN) 



2nX 2 N 2 

The asymptotic behaviour of Hankel function of complex argument leads to the following result 



B z = 



B z 



$ / nXN 
2irX 2 N 2 \^2r~ 



1/2 



for 



2irX 2 N 2 



In 



for 



r > A7V 



r < AiV . 



(78) 



(79) 



(80) 



It follows from (78) and (80) that in the centre of vortex the magnetic field tends to infinity. However, 
in reality the field can not be infinite and these formulae are not valid near normal core of vortex (r ~ £) . 
Then the expression (80) can be written as 



5(0) 



2^A 2 iV 2 



hiK . 



(81) 



6. Penetration of Magnetic Field into II Type Superconductor 

Assume that space-time admits a time-like Killing vector £ Q 

re* = -a , 



(82) 
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where —A = — iV 2 is the norm of time-like Killing vector, u a is the four-velocity of superconductor as a 
whole being parallel to the time-like Killing vector £ a 

u a u a = -1 , C = Nu a . (83) 

Then 4-momentum P a of the superconducting Cooper pairs is 

P a = mcu a s . (84) 

The four-velocity of superconducting electrons u™ can be decomposed in the form 

a — i — q .a 

(85) 



where vf is the relative velocity of the Cooper pairs. 
The energy of a superconducting pair is 

W = -P a C = mcu^ a = -mc " ?a = = N - . (86) 

y/1 - Vy C 2 

In the case of nonrelativistic motion, the kinetic energy of superconducting pair is 

2 

W = Nmc 2 + N^- . (87) 
The density of total kinetic energy of superconducting pair is 

W kin =N^L 2 . (88) 

where jf = n s evf is the superconducting current density 

Taking into account Maxwell equations written in vector form 

V x (NB) = -— j s (89) 
c 

one can get the expression for kinetic energy (88) in a more useful form 

W km = A^-(V x (NB)) 2 . (90) 
Sir 

According to expressions for the kinetic energy of supercurrent (88) and for energy of magnetic field the 
free energy of the whole superconductor is equal to 

^ sB = ^so + J [B 2 + A 2 (V x {NB)f]^dV . (91) 

According to Leibnitz formula 

(n^upA^) , a = (rj^upA^) B v + (if^ v u & B v , a ) A, , (92) 
Expressing vector potential A u through the magnetic field 

A„ = g^V ajTX u 7 B T , x (93) 

and using Stokes Theorem 13 

J (rT^upA^) , a ^—gdV = j {if^upA^B v ) n a dS = (94) 
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one could obtain the following expression for full energy of superconductor 



£ = ^- I B(B + A 2 V x Vx(NB))y/=gdV . (95) 



8tt 

Taking into account that magnetic field B is governed by Eq. (77) one can get 



£ = -tr m ■ (96) 



Inserting (81) into (96) one can get 

/ *„ \ 2 

N\nn . (97) 



AttXN 

For the superconductor embedded in applied external magnetic field the density of Gibbs free energy for 
the unit length of vortex 

g = S -N^. (98 ) 

takes the minimum value at the equilibrium. 

From this formula one can see that for the comparatively weak field Bq, Gibbs energy G > and the 
formation of vortex is impossible. However there exists the magnetic field B c \ from which G becomes negative 
and vortex formation becomes energetically favorable. It follows from (98) that 

The inner magnetic field should be found as a solution of the interior Maxwell equations. However for the 
superconductor embedded in an applied magnetic field B out , the interior magnetic field can be modelled as 
-So = B out N~ 2 for simplicity of calculations. Therefore, the free energy for the superconductor is 

F sB = F s0 + = F s0 + ^-N 5 . (100) 

When F s b = F n , the superconductor is transformed to the normal state and 

F n - F sQ = . (101) 

87T 

On the other hand according to Ginzburg-Landau theory 7 the free energy of superconductor in the 
gravitational field in the absence of external magnetic field is 

F s0 = F n + Na\iP\ 2 + N^\ 4 . (102) 

One can find the value of l^l 2 in which free energy has a minimum, by solving ^ffi = 0. Simple calculation 
gives 

|V| 2 = -a/(3 . (103) 
Inserting (103) to (102) and equalizing the result to (101) we will have 

9 Aira 2 

B ln = J^I ■ (104) 
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Then using the expressions for £, A and <&o we obtain for B, 

V2$p 
AttX£_N 

The second critical magnetic field 17 is 



Bom = 7^7772 • ( 105 ) 



B c2 = V2nB cm = . (106) 



7. Application of General Relativistic Electromagnetic Effects to Neutron Stars Interior 

The dependence of the coherence length and penetration depth on the gravitational field may be important 
for neutron star physics. There are two ways in which external or trapped magnetic flux can penetrate the 
proton superfluid of neutron stars. This depends upon two important lengths: the proton coherence length £ p 
and the London penetration depth A^. If £ p is larger than v^Al, there are small normal regions containing 
flux interspersed with field- free superconducting regions (type I superconductivity) . If instead £ p is smaller 
than ^/2\l, magnetic flux can penetrate the superconductor without destroying the superconducting state. 
Realistic estimations show that £ p is much less than the London penetration depth A^ for protons in the 
neutron star interior and hence the proton superconductor is expected to exhibit type II behavior. 

However after recent observation of long period precession of neutron stars it was pointed out in the 
literature, see for example Link, 14 that the axis of precession of some neutron stars may not be aligned with 
the axis of magnetic field. As the rotation of neutron stars causes a lattice of quantized vortices to form in the 
superfluid neutron state and it is generally believed that the proton superfluid is a type II superconductivity, 
which means that it supports a stable lattice of magnetic flux tubes in the presence of magnetic field, Link 
suggested that this two type of vortices might interact quite strongly due to the fact that the axis of the 
rotation and the axis of magnetic field are not aligned. 

Due to this reason it was discussed 18 that superconductivity inside neutron stars in fact may be type I 
because conventional picture of II type superconductivity follows from the standard analysis when only a 
single proton field is considered. In particular it has been shown the correlation length or coherence length 
£ should be replaced by the actual size of proton vortices 

5 = Z/y/i ■ (107) 

Since the parameter e arises from the strong interaction between proton and neutron superfluids, which 
is about 10~ 2 , the length (107) is much bigger than the usual proton coherence length. Consequently the 
Ginzburg-Landau parameter to be bigger than l/y/2, which in its turn causes that the superconductivity 
inside neutron stars to be type I. 

Naturally the question arises whether the general relativistic effects arising from the strong gravitational 
field of neutron star can change the type of superconductivity or not. As we have shown here, due to the 
gravitational effect the coherence length is also modified by factor of N (the value of the parameter N can 
reach ps 0.7 for the typical neutron stars), but it does not lead to change of the type of superconductivity, as 
gravitational field reduces the magnitude of the penetration depth by the same factor so then the Ginzburg- 
Landau parameter will remain unchanged. 

It is well known that in superconductors of II type below lower critical field B cl there will be complete 
expulsion of the field; above an upper critical field B c2 superconductivity will be destroyed; in the intermediate 
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range B c \ < B < B c2 the superconductor will allow the magnetic field to penetrate, not homogeneously but 
confined to quantized flux tubes. Each fluxoid carries a quantum of magnetic flux <!>o — hc/2e = 2 x 10~ 7 
Gauss • cm 2 . Within the core of a fluxoid (of radius £), the matter is in its "normal" state. The field strength 
rises to about B c i within this normal core. Around this the matter is in the superconducting state, and the 
field strength decreases exponentially away from the core with a scale length A^ . For protons in the interior 
of a neutron star, the coherence length is £ p = 0.6 x 10 _12 p* j^A" 1 cm and the London penetration depth 

is Al = 0.9 x 10~ n /9 p cm , where p Pj i3 is the mass density of protons in units of 10 13 g • cm~ 3 , and A p 
is the energy gap of the proton superconductor in MeV (Ref. 19). 

The vortex filaments consisting of neutron vortices and proton fluxoids inside superconducting region of 
neutron stars interact each other with an interaction energy which is equivalent to repulsion 20 



\4ttX 



2 

K, 



o © • (108) 



Due to the gravitational effect the interaction energy would be less when compared with what could be 
expected in the flat space-time by a factor 

n i n d N — l 

N- 7 / 2 e*— (109) 

which is tiny and may be neglected during accounting one of the main forces which has been acted on 
proton fluxoids. Here d = 5x 10 _10 (.Bi2) _1 is the spacing between proton fluxoids. But there are the most 
important forces that act on proton fluxoids such as magnus force, buoyancy force, drag force and tension of 
fluxoids. But there are the most important forces that act on proton fluxoids such as magnus force, buoyancy 
force, drag force and tension of fluxoids. The collective effect of these forces leads fluxoids to move radially 
outward. The force 21 

Jn — g , l liu J 

acting on unit length of fluxoid (due to neutron vortices) does not depend on gravitational field and remains 
unchanged. Here R c is the core radius, £l s and Q c are the rotation rate of the core superfluid and the rotation 
rate of the crust respectively. However the other forces such as buoyancy force arising from the existence 
of the magnetic stress in the core of a fluxoid and drug force that acted by the electron gas to a flux tube 
moving with the velocity v depend on the gravitational field and would be modified when compared with 
the Newtonian ones 22 

2 



= / 'I'd V I , f\ 



fb = — —In - , (HI) 

Jb \A-kXNJ R c U/ 

3tt n e e 2 $ 2 v p 

U ~ Q4E f (e)XN* c ' [ ' 

where Ef — NEf(e) is the electron Fermi energy in the stationary gravitational field. In spite of the fact the 
gravitational field of neutron stars does not change the value of these forces acted on fluxoids enough, the 
motion of fluxoids inside neutron star might be changed. Considering the equation of motion of fluxoids 21 
one can eventually see that difference of the rotation rate of the core superfluid from the rotation rate of the 
crust also depends on the value of the gravitational field. 
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Since the magnetic field decay from the core of neutron stars depends on interpinning of neutron superfluid 
and proton superconducting fluid 21 , then the slightly modified equation of motion leads to the change of the 
decay timescale, but this modification is very small, and could not be important for considering the lifetime 
of the magnetic field persisting in neutron stars. 

8. Conclusion 

Following up our previous research on the electromagnetic properties of superconductors in gravitational field 
we have considered some general-relativistic effects associated with the influence of stationary gravitational 
field on quasi-stationary electromagnetic effects in conductors and superconductors. 

By solving Maxwell equations for the time-dependent magnetic field penetrating inside the conductor 
we have shown that the skin depth in the conductors embedded in the gravitational field strongly depends 
on the gravitational field. However this effect is almost negligible for the laboratory conductors in the weak 
gravitational field of the Earth. 

The penetration depth for the magnetic field and the coherence length in superconductors also depend 
on the gravitational field and proportional to the lapse function of the gravitational object. However the 
parameter of Ginzburg-Landau remains unchanged and does not lead to the general relativistic modification 
of the type of the superconductivity. 

Two critical magnetic fields that destroy the superconducting state also depend on the gravitational 
field: the gravitational field increases the value of the critical magnetic field. Thus it has been shown that the 
strong gravitational field does effect to the quasi-stationary electromagnetic processes inside superconductors 
and one must take them into account during studying neutron star physics. 
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